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In higher Landau levels (N > 1), the ground state of the two-dimensional electron gas in a strong 
perpendicular magnetic field evolves from a Wigner crystal for small filling v of the partially filled 
Landau level, into a succession of bubble states with increasing number of guiding centers per bubble 
■ as v increases, to a modulated stripe state near v = 0.5. In this work, we compute the frequency- 

D ' dependent longitudinal conductivity a xx (u>) of the Wigner and bubble crystal states in the presence 

of disorder. We apply an elastic theory to the crystal states which is characterized by a shear 
^-j- ■ and a bulk modulus. We obtain both moduli from the microscopic time-dependent Hartree-Fock 

' approximation. We then use the replica and Gaussian variational methods to handle the effects of 

disorder. Within the semiclassical approximation we get the dynamical conductivity as well as the 
pinning frequency as functions of the Landau level filling factor and compare our results with recent 
microwave experiments. 

c3 

! PACS numbers: 73.43.-f, 73.20.Qt, 73.21.-b 
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INTRODUCTION 



In the presence of a strong quantizing magnetic field, the two-dimensional electron gas (2DEG) is expected to 
crystallize below a filling factor v ~ 1/6. 5i The resulting electron solid, or Wigner crystal (WC), has been intensively 
studied, theoretically, since the early days of the quantum Hall effect. On the experimental side, many groups have 
reported the observation of a strong increase in the diagonal resistivity p xx together with non-linear I-V characteristics 
and broadband noise when the filling factor of the lowest Landau level is decreased below 1/5. These observations have 
' been interpreted as the pinning and sliding of a Wigner crystal^ Early microwave absorption experiments 3 - detected 
l— ', a resonance in the real part of the longitudinal conductivity tensor, a xx (u>) , that was attributed to the formation of 
a pinned Wigner crystal. (A more recent experiment reports the observation of not one but two resonances at low 
filling factor. The origin of these resonances are not yet fully understood.) More recent studies show a resonance in 
\Q ■ the absorption at filling factors close to v — 1,2,3 where the formation of a Wigner solid is expected in very clean 
Q\ ' samples^ 

lO . In this paper, we concern ourselves with a series of microwave absorption experiments^ on the so-called bubble 
phases in Landau levels N — 2. The bubble phases, first introduced by Fogler and Koulakov^ are basically crystals 
made of clusters of M electrons. More precisely, the guiding centers of the cyclotron motion (the cyclotron radius 
is given by R c = \/2N + \i where I — y^Hc/eB is the magnetic length) of the M electrons are concentrated in a 
small area of radius Rm = V 2M£ at each lattice site. In the Hartree-Fock approximation, the 2DEG ground state, 
in Landau levels N > 1, evolves from a Wigner crystal (M = 1), near integer filling, into bubble crystals with 
M = 2, N + 1 as the filling factor increases, and terminates with the quantum Hall stripe state near half-filling of 
the partially filled leveliS. (The sequence repeats in inverse order for holes for partial filling factors between 0.5 and 
' ^ ! 1.0.) All these transitions are believed to be first order. 

Microwave absorption experiments* reveal a resonance in the second Landau level (N — 2) near integer filling, 
which was attributed to the formation of a Wigner crystal. The frequency of this resonance decreases, as the filling 
factor v — v to t — 2iV of the partially filled level increases, until it reaches a plateau at v — v\ « 0.16 where a second 
resonance appears (see Fig. 3 of Ref. 0). This second resonance is attributed to the formation of a bubble crystal with 
M = 2. The frequencies of both resonances are nearly constant between v\ and V2 ~ 0.26. At v-i, the first resonance 
disappears and the frequency of the second one starts to decrease. It has been suggested that the presence of the two 
resonances between v\ and V2 can be explained by the coexistence of the two crystal phases in this range of filling 
factor 

In a recent paper, we have studied the dynamics of the bubble crystals in the absence of disorder—. We have shown 
that these crystals have a low-energy magnetophonon mode with dispersion u> ~ fc 3 / 2 , at small wavevector, that is 
typical of a pure Wigner crystal in a strong magnetic field. In the presence of disorder, this low-energy magnetophonon 
mode becomes gapped, with lowest energy at the pinning frequency. In this paper, we compute this pinning frequency 
and the evolution of the longitudinal part of the dynamical conductivity, a xx (uS) , as the 2DEG makes transitions 
from the WC to the bubble states and compare our results with those of Lewis et ali (We have previously addressed 
the dynamics of the disordered stripe state in Refs. ^] and Il2t) . The pinning frequency is extracted from the real 
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part of the dynamical conductivity i.e. from -ft [a xx (q = 0, u)] computed in the presence of disorder. To deal with the 
disorder averaging, we use a combination of the replica trick and Gaussian variational methods (GVM)i" A similar 
technique has been used by Chitra et al^ to compute the dynamical conductivity tensor of the Wigner crystal in the 
lowest Landau level. While this method appears to give reasonable results for this and other pinned systems, some 
controversy has arisen regarding the width of the resulting resonant peak for a pinned Wigner crystal in a magnetic 
fieldi£*i&. We address this issue toward the end of this paper. 

To deal with the bubble crystals, some modifications of the replica method are needed. For example, electrons 
and bubbles are more extended objects in higher Landau levels and so it is necessary to include a form factor in 
the calculation of the dynamical matrix that enters the elastic action as well as the effective disorder potential. We 
find that the evolution of the bulk and shear moduli with filling factor is essentially driven by the changes in the 
form factor of the bubbles. A softening of the lattice appears at the critical filling v\ introduced above when the 
outer radii of adjacent bubbles touch one another and, as we will show, is correctly captured by a dynamical matrix 
computed with the correct form factor. Moreover, the gap between the two pinning resonances seen in the microwave 
experiment may be understood as being partially due to the change in form factor in the transition from the Wigner 
to the M — 2 bubble crystal. In our calculation, the bulk and shear moduli are extracted from the density-density 
response function computed in the Generalized- Random-Phase Approximation (GRPA)^ 

Our paper is organized in the following way. In Sec. II, we build an elastic model for the description of the pure 
crystal states and extract the dynamical matrix from the GRPA calculation. In Sec. Ill, we briefly review the replica 
and GVM that we use to handle the disorder. Our numerical results for the dynamical conductivity and pinning peak 
are presented in Sec. IV. We compare our results with the experimental data in Sec. V. Section VI is devoted to 
a discussion of the width of the resonance obtained in this method as compared to others, and we conclude with a 
summary in Sec. VII. A brief account of this work has appeared previously^ 



II. ELASTIC MODEL FOR THE BUBBLE CRYSTALS 



A. Form factor 



In describing the bubble crystals in higher Landau levels, we assume, as usual, that the electrons in the filled levels 
are inert. The electrons in the partially filled Landau level N are assumed to form a triangular crystal structure^ 
with M electrons associated with each lattice site. It is convenient to describe the crystal by the Fourier components 
(p (K)) of the guiding-center density where K is a reciprocal lattice vector of the triangular lattice. The guiding 
center density is related to the real electronic densitjiiS by the relation 



(n (K)) - J dre- lK r (n (r)) = N V F N (K) (p (K)) , 
where (n (K)) is a Fourier component of the real density, N v is the Landau-level degeneracy, and 

( L° N (x) is a generalized Laguerre polynomial) is the form factor of an electron in Landau level N. 
The Hartree-Fock energy per electron in the partially filled Landau level is given by 

l£[iT(K)(l-^o)-*(K)]Kp(K)>| 2 , 



(1) 
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where v — vtot — 2 N is the filling factor of the partially filled level. The Hartree and Fock interactions in Landau level 
N are given by 
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where k is the dielectric constant of the host material. Figure ^ shows the Hartree-Fock energy per particle for the 
bubble crystals in N = 2. The first order transition from the WC to the M = 2 bubble phase occurs at v\ = 0.22, that 
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FIG. 1: Energy per particle for the bubble and stripe states in Landau level N — 2. The region between the dashed lines is 
where the two resonances are observed to coexist in the microwave experiment. 



between the M = 2 and M = 3 bubble phases occurs at v% — 0.37 and that to the stripe phase occurs at ^3 = 0.43. 
In Fig. the energy is in units of e 2 /k£q where £q is the magnetic length evaluated at filling factor i^o = 0.1 using 

k£ V 2N + u \ K £ J ' [ ' 

an expression that takes into account the filled Landau levels. It is important to notice that the variation of the 
magnetic length £ through the M — 1 or M — 2 bubble phases is quite small for N = 2 so that we can effectively 
consider the variation of the filling factor to be equivalent to a variation of the density of electrons. In the microwave 
experiments^, two resonances are observed in the region bounded by the dashed lines in Fig. ^ which is argued to be 
the signature of the coexistence of the two phases. 

Once the (p (K)) s are known in the HFA, it is possible to compute the Matsubara two-particle Green's function 

X&,'& (k, t) = -N v (T T p(k + K, r) p(-k - K', 0)) , (7) 

in the GRPA 1 -. In Eq. (JJJ , p = p — (p) , k is a wavevector in the first Brillouin zone of the triangular lattice and %- 
is the time-ordering operator. The analytical continuation Xk'k' (k u °f ^ n ^ s Green's function gives the density 
response function whose poles are the collective (density) excitations of the crystal. The lowest-energy collective mode 
is the magnetophonon mode. We can get its dispersion relation, ujqrpa (k) by following the lowest-energy pole with 
non- vanishing weight of Xk.=o k'=o Ck u + ^) as ^ is varied in the Brillouin zone. 

In order to compute the dynamical conductivity <r xx (u>) of the disordered bubble crystals using the replica and 
GVM, it is convenient to construct an elastic theory valid for the low-energy excitations. As has been demonstrated in 
Refs. fllllTl we can do so by proceeding in the following way. We model the crystals as being made of rigid extended 
objects centered on each lattice site R with a density profile (or form factor) given by h (r). Because each bubble has 
M electrons, we have the normalisation 

J drh (r) = M. (8) 

The dynamics of these objects is described by the displacement fields u (R, t) at each lattice site R. The electronic 
density can therefore be given by 



n (r, t) = ^2 h (r - R - u (R, t)) 

R 



(9) 
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In the ground state, the HF energy given by Eq. (|3J is very well reproduced if we use the form factors 



h (q) = J dre-^h (r) 



-q 2 i 2 /2 



for the WC in TV = and by 

Mq) 



(l-(^) 2 + ±(^) 4 ), 



for M = 1 



(2 - \ [qlf) (l - (q£) 2 + § {qtf) , for M 



= 2 



(10) 



(11) 



for the bubble crystals in N = 2. These form factors are obtained from the Slater determinant for the wavefunction 
of a bubble of M electrons in Landau level N 



^ N (ri,r 2 , ...,r M ) 



¥n,i (ri) <^Ar,i (r 2 ) ... <pn,i (tm) 
</>iV,M-i (ri) <Pn,m-i (1*2) ■■■ <PN,M-1 ( r M") 



(12) 



where (r) = C w , m (r) 1 ™^ 1 e -7«V^-™) e 4™7^L 



is the normalized wave function of an elec- 

tron in the symmetric gauge A = (—Boy/2, Bqx/2, 0) with angular momentum m. The associated one-particle density 
is simply 



Mr) 



M 



.i=2 



II dr. 



-M-l 



(r,r 2 , ...,r M )| 2 = |^Ar, m (r)| z 



(13) 



In the presence of a transverse magnetic field B = —Boz, the equation of motion for the displacement field of the 
rigid object of charge Me and mass Mm* on site i is given by 



(Mm*) 



d 2 u (Ri,t) MeB du (Rj, i)xz 



d/ 2 



(14) 



where Z) (Rj — Rj ) is the dynamical matrix of the crystal (we use a "hat" to indicate that a quantity is a 2 x 2 
matrix). The dynamics in the strong-magnetic field limit is obtained by setting m* = in Eq. (|14fl so that, after 
Fourier transforming, we get 



where we have defined 



MHdu (k,i) 
~7 2 dt 



x z - D (k) • u (k, t) = 0, 



D(R i -R J -) = ^^B(k)e ik -( R -^ 



(15) 



(16) 



and 



u(k,i) = ^5>- ik ' R u(R,<) : 

VJVs R 



(17) 



with N s the number of lattice sites. Fourier transforming in time u (k, t) = J ^e luJt u (k, oj) , we have, in matrix form 



D xx (k) D xy (k) 



Dy X (k)+i^L0 



Dyy (k) 



Ua, (k,w) 
My (k, Uj) 



Using Eq. (|18|l . we get for the magnetophonon dispersion relation 



hM 



det 



D(k) 



(18) 



(19) 
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B. Euclidean Elastic Action 



Using Eq. ifTSjl. we can now write the Euclidean action So of the elastic model for the pure bubble crystal. We 
have 

S ° = ^fYl X! u a (k,uj n ) G^ 1 (k,uj n )up(~k,-uj n ) : (20) 

where the displacement fields obey the single Landau level dynamics 2 ^ 

[ Ux (R),u v (B!)]=if5 u , w . (21) 

From now on, we set ks = fi = 1. In Eq. H20|) . T is the temperature, u> n — 2nn/T is the bosonic Matsubara 

frequency, and the displacement Green's function G^l (k, iuj n ) — J^ T dTe luJnT (7~ T u a (k, t) up {— k, 0)) So is related to 
the dynamical matrix by 

(kiu,)- £4 ( , r D w ( k ) ^-D xy (k)\ 

( ' lUJn) - AP^l+ul) Dy.<k) D xx (k) )■ (22) 

Once the Green's function has been determined in the presence of disorder, we can easily compute the conductivity. 
Since the electric current is carried by the charge, the current density can be expressed as: 

.Me du (k, r) 

j(k ' T) = l ^7^^- (23) 



The conductivity is then determined by the Kubo formula to be 

(24) 



l/T 

dre^ T (j a (k = 0,r)j> (k = 0,0)> 



M 2 e 2 

i W G r 4(k = 0, W ), 



where 5 is the area of the 2DEG and v c is the unit cell of the bubble crystal. Since D (k = 0) = 0, it is easy to check 
using Eqs. i|22|) and (|24|l that, in the pure limit and in the strong magnetic field approximation, the electromagnetic 
response of the system is purely transverse for all frequencies oj « u> c ; i.e., we have in this limit 

/ nec \ 

?H = (1 g J- ( 25 ) 
where n = M/v c is the electronic density. When the mass term is not neglected in Eq. (|14H . we have instead 

?M= ^^f ( iW "M, (26) 

and the absorption (which is proportional to [a xx (to)]) is at the cyclotron frequency lo c as it should be from Kohn's 
theorem 2 ^. 



C. Dynamical matrix in the GRPA 



In this work, we obtain the dynamical matrix from the GRPA. This procedures, which uses the density response 
function XkK' 0*' w ^) f° un d microscopically to obtain the dynamical matrix for the elastic theory, was developed 
to study pinned quantum Hall stripes, and is described in detail in Ref. Il2t Briely, we use the fact that the density 
fluctuations are related to the displacement field via the relation 

Sn(k + K, t) w -i^/N~ s h(k + K) (k + K) • u(k), (27) 

to obtain the relation 



(28) 
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between the bare displacement Green's function G^l (k, u) and the density reponse function. The density response 
function in the GRPA depends on the guiding center densities (p (K)) s computed in the Hartree-Fock approximation. 
The components of the matrix (k, r) (or of Dgrpa (k)) are obtained from both the poles and the weight of the 
density response function. 

It is interesting to compare the dynamical matrix obtained in the GRPA from that obtained with other approxi- 
mations used in the literature. If the filling factor v is not too small, a rapidly convergent expression for the classical 
dynamical matrix of a lattice of rigid objects with form factor h (r) is given by 



-^-'classical t^KJ — 



K 



^ k + K )' 2 (k + K)(k + K)-^KK 



(29) 



|k + K| ^ |K| 

where n s = N s /S is the density of crystal sites. This expression is not suitable, however, for a lattice of point particles 
k 2 . 



(p.p.). In this case, we can use the expression of D pp (k) given by Bonsall and Maradudm 20 which is valid to order 



To compare the three expressions for D (k) , we consider the parameterization of the dynamical matrix in terms 
of the Lame coefficients A and p. The deformation energy of the crystal can be expressed in terms of the dynamical 
matrix as F = 5 ^2 a g u a (k) D a g (k) up (— k) . On the other hand, the same energy is written^ in elastic theory 

(and for a triangular lattice) as F = | J dr A (e^) 2 + 2pe 2 ^ , where e a ,g — ^ ( Hgf^ + a ^p( r ) j } s t ne strain tensor. 

The coefficient p is actually the shear modulus of the lattice while the bulk modulus is given by B — = A + p. 

Comparing the two expressions for F, it is easy to see that the dynamical matrix, in the long-wavelength limit, can 
be written in terms of the Lame coefficients as 

F><x,p (k) = n^ 1 (A + p) k a kf3 + n^pkfSa^. (30) 

The long-range Coulomb force, absent in the elastic theory, is very important for correctly capturing the low-energy 
physics. Its inclusion leads to an additional term in A, i.e. to 

A ^HH-+> ? - (3D 
The long- wavelength dispersion relation of the magnetophonon mode is given (to order k 2 ) by 

' e2 -^/(\ + 2p)pk 2 . (32) 



fiM 



It depends mostly on the shear modulus p and on the long-range Coulomb part of A. We will show later that the 
pinning behavior is also mostly determined by these two quantities. 

Using Eq. I|3U|I , we can obtain pgrpa from our numerical data for Dgrpa (k) . The shear modulus /^classical can be 
obtained by making a small wavevector expansion of Eq. (|29|l taking into account the appropriate form factor h (q) 
that depends on M and on the Landau level N. To order fc 3 , we find that the Lame coefficients can be written as 

/27m 2 e 2 \ (M 2 a 4 \ 

Aclassical = I — - — I I + «o + + 2a 2 - a 3 + — + (3 k J , (33) 

/ 27rn 2 e 2 \ / a 4 \ /0 

^classical = I ~ I ^3 + ~^ j ■ (34) 

The detailed forms of the coefficients a^s and /3^s are given in Appendix A. 
For a triangular lattice of point particles, we have instead^ (to order k 2 ) 

Vp.p = ~n a (« + 2b) , (35) 
Pp.p = n s b, (36) 

where a = -1.225323^^,6 = 0.245065-^si, with v c = ^-al the unit cell area of the bubble crystal with M 
electrons. 

In Figs. [3 and we compare the three different approximations for the shear modulus. As can be seen from Fig. 
13 a point lattice is a good approximation in Landau level N = (in comparison with the more exact GRPA result) 
because the shear modulus does not vary much with filling factor. In JV = 2, however, the form factor must be 
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FIG. 2: Shear modulus computed in different approximations for Landau level N = 0: for a classical lattice with form factor 
h(G) (full line); for a lattice of point particle (dashed line) and in the GRPA (line with squares). 



0.007 



0.006 - 



0.005 - 

CO qW.wv^ 

y 

CM 

a> 0.004 - 
o 
<n 

■- 0.003 - 



3_ 0.002 - 



0.001 - 



0.000 



b Classical with form factor 
■ classical {point particles) 
— GRPA 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 

V 



FIG. 3: Shear modulus computed in different approximations for Landau level N = 2: for a classical lattice with form factor 
h(G) (full line); for a lattice of point particle (dashed line) and in the GRPA (line with squares). 



considered to account for the softening of the lattice near the transition between two crystal phases. In this case, 
the GRPA result is qualitatively similar to that of the classical approximation although, quantitatively, the quantum 
corrections are important especially when considering the jump in the shear modulus at the transition. In the small 
v limit, all three approximations give the same result for /i as expected. 

Note that in Fig. the shear modulus decreases near 1/2 where the Hartree-Fock calculation predicts a transition 
to another bubble phase with M = 3. In the transport experiments of Lilly et ali^, there is no sign of such transition 
and, in the rest of this paper, we will refer to the transition region as the region around filling factor v\. The range 
of filling factor where the M = 3 is stable is quite small in the HFA and, in a more accurate calculation, this phase is 
probably unstable with respect to the stripe phase. Indeed, the M = 3 bubble state is absent in DMRG 25 calculations. 
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III. REPLICAS AND THE GAUSSIAN VARIATIONAL METHOD 

A. Disorder Coupling 

Now we take into account the disorder effect. We assume that the main source of disorder comes from the interface 
roughness and that it can be represented by a spatially uncorrelated Gaussian random potential V(r). The disorder 
action reads 

SU P = J dr dTV(r)n(r,T), (37) 

where V (r) has a Gaussian distribution function which leads to the correlator 

V(r 1 )V{r 2 )=V<?v c 5(r 1 -r 2 ). (38) 

Here the overline denotes an average over disorder configurations. 

The electron density operator n (r, r) in Eq. 1)3 7|l must be approximated more accurately than was needed in in 
order to capture the possibility of pinning by disorder. Following Giamarchi and Le DoussaU^, under the assumption 
of small Vu(r) (which is justified for weak disorder) we write for the bubble crystals 

n (r, t) = n - n V • u (r, r) + n s ^ h(K)e lK ^- u( - r - T ^ , (39) 
where n a = n s h (0) is the average electronic density. Fourier transforming, we get 



n (q, t) 



iV <S qi o - i y/N~ s Mq • u (k, t) + n s ^ h(K) f dr e «-[r-u(r,T)]-<q.r 



(40) 



where q = k + K. Only the last term in Eq. 140(1 captures the short wavelength oscillations in the charge density 
and allows pinning by impurities^ so that we will drop the first two terms of this equation when using n(q, r) in the 
disorder action. 

B. Effective Form Factor for the Coupling with Disorder 

To simplify the coupling with disorder, it is usual to keep only the summation over the first shell of reciprocal 
lattice vectors in Eq. (|4t)(l (i.e. 6 vectors for a triangular lattice) as these give the strongest restoring force for the 
displacement fields. Using Eqs. Q and 0, one would then be tempted to use, for the form factor h(K) that enters 
Eq. (|40[l. the expression 

h(K) = ™F N (K)(p(K)), (41) 

where (p (K)) is the guiding-center density computed in the Hartree-Fock approximation. This choice, however, 
underestimates the fluctuations in the potential as an electron wavepacket or a bubble moves across the disorder 
potential. We can capture this energy scale with an effective form factor, that we call H(K), in the following way. 
The density in the pure crystal is exactly given by 



(r) = l^n(K) e lK ' r . (42) 



S K 



The fluctuation in the disorder energy when a large patch of the crystal is slid uniformly is then 

(U 2 ) = j dr I dr'Sn (r) Sn (r') (v (r) V (r0\ , (43) 
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Using Eqs. ©, {TUJ , and JTJ we find 

TT 2 \ = l 

s 2 



(U 2 ) = ^ElWI 2 ' (44) 



K/0 



mp — 



= ^^El^( K )l 2 KMK))| 2 , 

where (p (K)) is computed in the Hartree-Fock approximation. On the other hand, if we describe the fluctuation in 
the disorder energy by Eq. I|4l)|l . keeping only one shell of reciprocal lattice vectors with modulus Kq, we get for the 
same energy 

AT V 2 

(U 2 )*^6\H(K )\ 2 . (45) 

We define our form factor, H (Kq) , by requiring both expressions to be equal so that 

1 M 2 

\H(K )\ 2 = E \ F " ( K )! 2 K<°( K )>| 2 ■ ( 46 ) 

In practice, we keep up to 50 shells of reciprocal lattice vectors in the summation on the right hand side of Eq. i|46[) . 

From now on X^k^o m ^ ne disorder action will mean a summation over the first shell of reciprocal lattice vectors. 
This simplification allows us to compute the Green's function in the presence of disorder in a relatively straightforward 
manner while retaining the essential physics of pinning so that our results are qualitatively correct. The major effect 
of these approximations is to replace the soft cutoff in wavevector that would enter through the form factor with a 
hard one in the reciprocal lattice sum. With this approximation, the impurity action with which we now work is 

f drdTV(r) \H (Kq) | ^ e iK-[r-u(r,r)]_ (47a) 

C. Replicas and the GVM 

A detailed account of the replica and GVM can be found in Ref. In previous work, we generalized this method 
to study pinning of quantum Hall stripes near half filling in higher Landau levelsii*^. Since the formalism can be 
applied with minor changes to the pinning of the bubble crystals, we refer the reader to these papers for details. In 
this section we only outline the procedure to obtain the so-called saddle point equations (SPE's). 

With the replica trick^i, one creates n copies of the original action, computes the replicated partition function Z n , 
and performs the disorder average on Z n . Also by taking into account the effective form factor as we discussed before, 
the resulting effective action is 

S eS = S^ + S^ p \ (48) 
1 n 

5 o Cff) = ^EE E «S(k,w n )G(g- 1 (k sWw )uj|(-k,- WB ), (49) 

a—1 k,^™ a,0—x,y 

sS - -<w E / dr i / dT2 E E \ H (^)i 2 ( 5 °) 

a,b=l •'° ^° r K^O 

xcos [K • [u a (r, ti) — u b (r, r 2 )]] , 

where v- lmp = V 2 v 2 n 2 , and a, b are replica indices that run from 1 to n. In obtaining the last line of Eq. I|50() we have 
neglected some rapidly oscillating terms. 

In the pure limit, the action is diagonal in the replica indices. Disorder averaging introduces a coupling among the 
replicas through the impurity coupling S^J in Eq. (|50|l . This coupling is non-Gaussian so we next apply the GVM 
to get a simplified expression. The fundamental idea of the GVM is to replace S c g with a variational action 5 var that 
is quadratic. We write 
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with the coefficients (G (k, u> n ) chosen to best match the original problem. This is accomplished by minimizing 

a free energy^ F var = F + T [(S) s — (<S , va r)g ] where F is the free energy associated with the action S val , and 
(■ • • )g nr indicates a functional integral over displacements, with S val as a weighting function. In Eq. (|51|l . G^b (k, uj n ) 
is the displacement Green's function, 

G$(k,u,„)= / 1/T dr(T r <(k,r)^(-k,0)) s . (52) 
Jo 

It is convenient to write it in terms of the bare Green's function via the equation 

(G- 1 )^ (k,c„) = G^" 1 (k,w n ) <5 afc - C$ k) , (53) 

where C^jg (ui n ) is the element of the variational self-energy matrix £. Note that there is no k dependence in £ because 
we have chosen our impurity action to be local in space; this will become clear when we find the SPE's below. Note 
also the obvious symmetries G ab = G ba and ( ab — ( ba . 

D. Saddle Point Equations (SPE's) 

Finding the extremum of the free energy leads to the saddle point equations (SPE's) 



QK) = 4« imp J^' T drUl- cos w n r) V' aB [B aa (r)} + £ [B ab (r)] 1 , (54) 

Cl Ma) K) = -4 Uimp f 1/T dr wsuwrV^ [B ab (r)} , 
Jo 



where 



and 



<f> [B ab (r)] = W (K )\ 2 K a K f3 exp 



1 



2 ^ K^B^r) 



(55) 



(56) 



S^(r) = <T r K(r,r)-^(r,0)] 2 ) Svar (57) 
= T ^ E K(k^„) + G bb p (k,Lu n ) - 2cos(w n r)G$,(k, Wn )] ■ 

It is apparent at this point that the self-energy has no k dependence. Moreover, if we assume that reflection 
symmetry for the bubble crystals are not spontaneously broken after disorder averaging, then the solutions of interest 
to Eqs. (|54|l and l|55|) will satisfy Q.y = and Q b = Q b . Using the symmetry [see Eq. (|30H ] of the dynamical matrix 
D xx (k x ,k y ) — D yy (ky,k x ) for a triangular lattice, we can also show that (® b = C,y b is a solution of the SPE's. We 
need thus solve for one self-energy component only. 

At this point, we take the limit n — ► in the replica index. In taking this limit, one assumes that the self-energy and 
Green's function matrices can be written in a "hierarchical form"—. A consequence of this is that the replica index 
becomes a continuous variable u € [0, 1] (The variable u must not be confused with the displacement held here.) For 
example, the self-energy matrix is now characterized by a diagonal component, £ Q , and an off-diagonal component, 
( a (u), such that 

(aa — * Cq, (58) 

C^ a) -Ca(«). for 0<«<1. (59) 

Similarly, G^ — > G a p, — > G a p(u) (0 < u < 1). Since the disorder potential F(r) is time independent, a 
further simplification one hnds is that the off-diagonal replica components Cck*^ an d £7°^ °^ are r independent^ 
so that G(k, u> n ,u) and £(k, U! n ,u) are different from zero only for lu„ — 0. We have 

G(k,u n ,u) = G{k,u)5 Un , , (60) 
C(w n ,u) = C(«)<W,o- ( 61 ) 
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The SPE's H54fl and (|55fl for the diagonal component Cxx — Cyy = C may now be written as 

f 1 f 1/T r~ 

CK)= / du({u)+4v imp dr (1 - cos (w„r)) V' xx B{t) 



where, from Eq. (|57|) . 



C(«) 



T 



S xx (r) = 2T ^ (1 - cos(w n T))G xa: (k,w n ), 



G xx (k,u) 



(62) 
(63) 

(64) 
(65) 



Note that Eq. (|rJ2)> also gives us C(w„ = 0) = f Q du((u). 

To compute the dynamical conductivities in Eq. <|24|l . we first need to solve the SPE's for the finite- frequency self 
energy C( w n)- If we take the analytical continuation of Eq. (|62l) . we get 



C rct H = e - 24t> imp \H (K Q )\ 2 K 2 I dt{e wt - 1) 



I K 2 [°° 
x Im exp < - / d/Im 

n In 



L k 



(1 - e -tf*) 



(66) 



where eo = C rct (0 + ) is a constant which will be determined below. Eqs. 1(60(1 and 1(61(1 indicate that for non-zero 
frequency u> n , 



G xx (k,u n ^ 0) = [G^-^k,^) - CK) 



(67) 



Thus, once we have obtained the self-energy, we can easily compute the finite-frequency conductivities in Eq. 1(24(1 by 
analytically continuing Eq. I|67|l to real frequency. The retarded Green's function for the displacement field is given 
by ' 



G^(k, o, ^ 0) = G™-\K u) - C rct (w) 



(68) 



and so 



G(k,w) 



D ra (k) - C rot (w) 



'2/2/ 



D xy (k) 



Du* (k) A™ (k) - C rct (w) 



det [L>] - C rct (w) tr [L>] + ( C rct (w„) ) - hPuj 2 /^ 



(69) 



Inserting this result into Eq. (|24|l . we obtain our final expression for the longitudinal conductivity in the presence of 
disorder, 



o~ xx {w) = CTyy(L0) 



C ct (uj) -M 2 cu 2 /£ 4 



(70) 



E. Semiclassical approximation 



Equation 1)66(1 for the self-energy is a very complex self-consistent equation. In studying the stripe phas o 11 ' 12 , where 
a depinning transition may take place and relative fluctuations in the displacement field can become very large, we 
found out that it was necessary to solve this equation exactly. In the case of the bubble crystals, the transition 
between the bubble phases are first order and we expect that the fluctuations in the relative displacement of the 
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bubbles are small enough so that we can use a semiclassical version of Eq. I|66(l . The semiclassical expression for the 
self-energy^ is obtained by expanding the exponential in Eq. I|56f) and keeping only the leading term. This leads to 



C rct H = eo + ^E feCk. " ^(k, u, - 0+)] , (71) 
s k 

where 

A = 12 Uimp |ff(if )| 2 ^o- (72) 

The semiclassical approximation is a powerful simplification when it is valid. In particular one sees that Eq. I|71|) 
is local in the frequency lo so that £ ret (o;) can be determined one frequency at a time. As we will see later, eo ^ is 
an energy offset that, in a pinned state, opens a gap in the phonon spectrum. 



A Constraint for eo 



The SPE H71(l becomes a set of two coupled equations for the real and imaginary parts of the self-energy if eo 
is known. Formally, eo needs to be determined self-consistently by solving Eq. (|71|) together with those for ((u). 
The SPE's have a family of solutions parameterized by u c where u c is, in principle, determined by minimizing the 
free energy. In the case of spatial dimension d > 2, u c determined in this way leads to 3? [tx QCt (w)] ~ lo 2 at small 
lo. This behavior is consistent with arguments by Mott as well as with some exact solutions^ (up to a logarithmic 
correction). For d < 2, however, this way of finding u c can yield an unphysical result in which, in the pinned state, 
the conductivity shows a true gap i.e. 3? [a (to)] vanishes below some finite frequency. For d < 2, what can instead be 



done is to impose the condition 3i [cr aa (a;)] 



at small us. This common procedure, although not fully understood, 



leads to physically reasonable results and we will adopt it in the computations that follow. (The reader should keep 
in mind that it leads to a broader response in the WC regime than is found by other methods. We discuss this issue 



C rct H 



lo at small lo and this guarantees 



below.) From Eq. (|70|l . this condition is equivalent to the condition S 

that the magnetophonon mode density of states vanishes at zero frequency as it should for a pinned system. 
To determine eo, we thus write, for small lo, 

(ret M ~ 6q+ IpLO. 



(73) 



With this constraint in Eq. I|71|l . using Eq. (|69|l and the symmetry relation ( x c x (lo) = Q C y(co), we can easily show that 
the condition for nonvanishing (3 leads to 



det [D] + e tr [D] -e\-\ (tr [D]f 
(det [D] - e tr [D] + eg) 2 



(74) 



Eq. 174|i is the constraint which we solve numerically to get eo- 



IV. NUMERICAL RESULTS 



We now discuss our numerical results for the behavior of the dynamical conductivity a xx (lo) with filling factor in 
Landau level N = 2. Figures 01 and [S] show the real part of the dynamical conductivity 3? [a xx (lo)] for different values 
of the filling factor v in the M = 1 and M = 2 bubble phases. At the transition between two bubble phases, the 
density is such that the cyclotron radius R c = \/2N+l£ of adjacent bubbles touches one another. Away from the 
transition region, however, the size of the bubbles is smaller than the lattice spacing. For all filling factors here we 
assume the bubbles are much bigger than the correlation length of the disorder (which we take as zero) . 

We plot in Fig. (a) the peak (or pinning) frequency f p k — 1-klo^ , (b) the full width at half maximum of the 
resonance, Af p k, and (c) the quality factor Q — f p k/ Af p k- Although we do not expect our theory to be quantitatively 
accurate, we choose the disorder level in our calculation so that our pinning frequency at v = 0.1 is equal to that 
measured experimentally^ at the same filling factor; i.e., f p k ~ 1.5 GHz with an electronic density of n = 3.0 x 10 11 
cm~ 2 . We find that our results are in a very good qualitative agreement with those of Lewis et ali£ with the exception 
of the region between v a = 0.16 to Vf, = 0.28 where the two phases are assumed to coexist. 

The variation of the conductivity with v is the same for the M — 1 and M = 2 phases. In both cases the conductivity 
curve has a non-Lorentzian lineshape with a peak frequency and width that decrease with increasing filling factor v (or 
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FIG. 4: Dynamical conductivity for the Wigner crystal in Landau level N — 2. 
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FIG. 5: Dynamical conductivity for the M — 2 bubble crystal in Landau level N = 2. 



density) and a peak height that increases with increasing filling factor. This behavior is qualitatively consistent with 
that of a weak-pinning model. Quantitatively, however, our pinning frequency decreases faster (we find f P k ~ l/^ 3 , 
for M = 1) than what is seen experimentally where the dependence f p k ~ l/^ 3 / 2 is found. The width of the pinning 
peak and the quality factor that we get, however, are typical of what are measured experimentally. 

If we compare the behavior of the pinning frequency in Fig. [5] with that of the shear modulus calculated in Fig. [3] 
we see that the pinning frequency behavior is opposite to that of the shear modulus within the whole range of filling 
factor considered. A more rigid crystal is less sensitive to the disorder potential and has, as a consequence, a smaller 
pinning frequency. This is precisely what is expected from collective pinning theorySl. 

The increase in the pinning frequency that we get after v w 0.18 in the M = 1 phase is not seen experimentally. 
Instead, two pinning frequencies separated by a gap of the order of 0.2 GHz are detected in the region between 
the dashed lines in Fig. ^ These two frequencies are found to be nearly frequency independent in this range {i.e. 
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from v a = 0.16 to = 0.28) where the two crystal phases probably coexist. Interestingly, the difference in pinning 
frequencies at the endpoints of this interval, f p k (y a ) — f p k (fb) > is 0.42 GHz in our calculation. Wc remark that the 
choice of the form factor H (Ko) [see Eq. i|4fi[l] that enters in the coupling with the disorder [Eq. 1)50(1 ] is important 
in getting the right order of magnitude for the drop in the pinning frequency. We find, for example, that replacing 
the X^k^o 1-^ C^o)| (• ■ •) in Eq. (|50() (the sum is here on the first shell of reciprocal lattice vectors only) by the more 
naive form X^k^o IMK)| 2 (• ■ ■) (where the sum is over all reciprocal lattice vectors) and using a cutoff in K in this 
last expression overestimate the sudden reduction in the pinning frequency seen experimentally^. The formalism we 
use needs to be generalized in order for us to properly compute the dynamical conductivity in the coexistence region. 

In Landau level N = 0, our numerical calculations give results similar to those of Chitra and coworkers^. For 
uncorrelated disorder, the pinning frequency increases continuously with decreasing filling factor (at fixed density). 
For a disorder level {vi mp in our calculation) corresponding to the typical peak frequencies observed experimentally 3 , 
we get a quality factor Q of the order of l.We have found that if one solves the GVM equations for very small disorder 
values in Landau level N — 0, a sharp resonance emerges with Q ~ A -1 / 4 . However, this behavior appears only at 
disorder strengths yielding resonance frequencies far below those observed in experiment. 

We note that experiments on higher density hole systems 3 ^^ 1 " have revealed dramatically narrower resonances, a 
phenomenon that our present calculations do not seem to reproduce. We discuss this issue in the next section. 



V. LONG RANGE COULOMB INTERACTION AND WIDTH OF THE PINNING PEAK 



One issue regarding the lineshape of the dynamical conductivity near the pinning frequency involves the observation 
in hole systems, in the Wigner crystal regime, of pinning peaks with quality factor Q of order 5 and even higher 30 31 . 

The possibility of such dramatic narrowing was suggested by one of usi£ as being a result of the long-range Coulomb 
interaction, in conjunction with the effect of the magnetic field. In Ref. fTR a simulation of an electron system in a 
strong magnetic field with a simplified disorder model yielded a remarkably sharp resonance, with width so narrow it 
could not be determined numerically. That this effect was related to the long-range nature of the Coulomb interaction 
was demonstrated by running the same simulation with a screened interaction, which yielded a more typical pinning 
peak with Q ~ 1. 

How such a narrow resonance might be attained can be understood in a self-consistent Born approximation (SCBA). 
Writing the self-energy as C,(oj) = z\ + iz2, the imaginary part Zq,(uS) at the resonance frequency sets the width of the 
pinning peak. The SCBA can be written in the form 

C( W )=A|A^(q, U ), (75) 

where A is the scale of the disorder interaction defined in Eq. (I72|) . The similarity of this equation with the 
semiclassical approximation in the GVM [Eq. 1|71|) ] is apparent, and indeed if one takes the imaginary part of Eq. 
(I75|l it is identical to the imaginary part of Eq. l)71f) . This equation can be written explicitly in the form 



j2 



,, \ (tr [D] f - dot [D] W - zjtr [D] + zj + z\ 

q [det [D] +z\ -zI~uj 2 - zxtr [£>]] 2 + z\ [tr [D] - 2zi] 2 ' 



A solution to Eq. (|76[) appears at first challenging for small A because the right hand side appears to vanish as 
A — > whereas the left hand side remains at one. A solution is always possible for any given value of z\ because the 
integral is logarithmically divergent as Z2 — > 0. This can be seen by noting that at small q, the trace and determinant 
respectively have the forms tr [D (q)] = Yq and det [D (q)] = aq 3 . (Note for short range interactions the powers of q 
in each of these expressions increase by 1.) If we choose lo — —zi(u>), i.e., the resonance frequency, the small frequency 
limit of the integrand gives 

^W^WUr (77) 



with q c = y z\Y/a. We then arrive at a solution of the form 

Z2 „ ^ e -^A (78) 

A very similar result was obtained in Ref. using a different method. The exponential behavior in Eq. I|78|) allows 
for a very narrow resonance, as was found in the simulation. As discussed in Ref. Il5l the result may be understood 
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FIG. 6: Numerical results for the bubble phases in Landau level N — 2 for (a) the pinning peak; (b) the width of the pinning 
peak; (c) the quality factor. 



as being due to a suppression of the collective mode density of states at low frequencies, arising from the effectively 
very large stiffness that is a result of the Coulomb interaction. 

Given the similarity between Eqs. I|76|) and H71jl. it is surprising that the GVM approach does not obtain this 
narrow resonance. The reason for this difference can be understood by the way in which the GVM is different from 
the SCBA, specifically by how eo is determined. By definition, eo = z\(lo — > 0), which is obtained in the SCBA by self- 
consistently solving the real part of Eq. I|75|l . A major limitation of this approach is that it does not correctly capture 
under most circumstances the competition between elasticity and pinning energy that sets the pinning frequency 26,27 . 
In the GVM, eg is obtained by solving Eq. J73}, which followed from requiring 3? [cr QQ (o;)] ~ uj 2 at small w. There is 
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a very close connection between Eqs. (|74|) and (|77|l : the integrands are identical for q/q c >> 1. Because q c becomes 
very small for small disorder, by noting that zi(u)) is not very different than eo near the resonance we can see the 
choice of eo guarantees that Eq. 1)76(1 will very nearly be satisfied from the portion of the integral where |q| > q c .The 
small differences in the two equations can then be accomodated with a value of Z2 not very different than zi.This 
leads to the more typical result for a pinned system, Q ~ 1. 

We conclude this section with a few comments. First, it should be noted that the very sharp resonances predicted 
by the SCBA and related worlsi£i±& are not what is observed experimentally. While the hole experiments yield peaks 
with a surprisingly large Q, they are many orders of magnitude smaller than what is predicted by Eq. (|78|) . If 
such a very narrow resonance is the correct expectation for the pinned Coulomb system, then there is presumably 
some dissipation mechanism that has not been identified. One cannot help but notice, however, that the results of 
replicas+GVM method are far more similar to the experimental results than the SCBA by itself gives. And while 
the simulation reported in Ref. suggests a very narrow resonance, it is possible that the simplified disorder model 
used there may give a different result than a generic disorder model (although previous studies have suggested this 
should not be the case^.) 

Secondly, it should be emphasized that the choice of eo in Eq. 1)74(1 really is very special. If eo is allowed to vary 
even by a very small amount from the solution of this equation, we have verified that the semiclassical approximation 
to the SPE's [Eq. 1(71(1 ] results in a very sharp resonance as expected from the SCBA. This demonstrates that the 
important difference between the method used in this paper and the SCBA is not really in the use of replicas or the 
GVM, but rather the introduction of the extra constraint. Indeed, it was suggested by Fogler and Hus©i& that if one 
can instead impose a condition 5ft [cr Q a(^)] ~ u x with a large value of x, then the Q of the resulting resonance is again 
extremely large - smaller than what is expected from the SCBA, but still orders of magnitude larger than what is 
seen in experiment. 

A third important point is that, although these considerations leave us with questions about the Q of the resonance, 
it seems likely the resonance frequency itself is well-represented by the calculations described in this work. This is 
because the pinning resonance is determined by collective pinning, which balances the lowest two energy scales - 
elastic and disorder interactions - to determine the restoring force for motion of the system as a whole. (Again, we 
emphasize that even a very small adjustment of eo results in a narrow resonance. The resonance frequency does not 
need to change noticeably to dramatically increase Q.) A resonance frequency determined by collective pinning is 
consistent with what was found in the simulation of Ref. 

As we have seen, the GVM+rcplicas method, together with the eo constraint, yields results that generally compare 
favorably with experiment. Other methods such as the SCBA yield much sharper resonances than are observed. 
Nevertheless, it is possible that such sharp resonances are the correct response for the model we are analyzing, in 
which case the mechanism that leads to pinning peaks with Q ~ 1 in experiment remains unknown. We leave these 
questions for future research. 

VI. CONCLUSION 

We have studied the behavior of the pinning peak of the Wigner and bubble crystals in Landau level N = 2 as a 
function of the filling factor. We used an elastic action obtained from the Hartree-Fock and time-dependent Hartree- 
Fock approximations and treated the disorder averaging using the replica trick and Gaussian Variational Method. 
Comparisons with recent microwave experiments in N = 2 show that the predictions of our model for the dependence 
of the pinning frequency and width of the pinning peak with filling factor compare favorably with experiment. At the 
moment, however, our approach can not reproduce the pinning frequencies observed in some range of filling factor 
in N — 2 where the Wigner crystal and the M = 2 bubble phases are believed to coexist. Finally, we discussed the 
relation of this method with those that predict a very narrow resonance. 
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APPENDIX A: LONG WAVELENGTH EXPANSION OF THE DYNAMICAL MATRIX 



The dynamical matrix is given by Eq. I|29() i.e. 

2wn s e 2 

K 

where 



D (k) = V [T (k + K) |A (k + K)| 2 (k + K) (k + K) - T (K) |A (K)| 2 Kk| , (Al) 



A(k) EE /l 2 (k), 



T(k) 



1 



|k + K| 



(A2) 
(A3) 



We want to find an expression for the dynamical matrix in the long- wavelength limit valid up to order k 3 . For the 
form factor, we use Eqs. (|1L)|) . (|1 1| > . We expand the interaction and the form factors to order A; 3 and define in this way 
the coefficients To, Ti, T2, T3 and Aq, Ai, Aj, A3 : 



T (|k + K|) » T (K) + Tj (K) (k • K) + T 2 (K) k 2 + T 3 (AT) (k • K) 2 
A (|k + K|) « A (K) + Ai [K) (k ■ K) + A 2 (K) k 2 + A 3 (K) (k • K) 2 



(A4) 
(A5) 



In these expressions, all coefficients depend only on the modulus of K. For K = 0, we need the expansion to order k 2 
of the form factor so that we define the coefficients Ai and A2 by 



A(k) = M 2 + Aifc + A 2 fc 2 



(A6) 



The coefficients Ai,A 2 and Ao,Ai,A 2 ,A3 depends on the particular form factor that is choosen while To (AT) 
Using the inversion and reflexion symmetries of the triangular lattice, we find, after some algebra, that 



D xx (k) 

Dyy (k) 
D X y (k) 



27m, e 2 



A£2 
M 2 



+ a + a.\ + 2a 2 + a 3 + CK4 + Pok J k 2 + (a 3 + ^ j 



Ky 



+ a + a>i + 2a 2 + a 3 + a 4 + Pok J k 2 + (a 3 + k 2 



2im s e 2 \ ( M 2 2a 4 n , , 

— - — I I —j— + u + ai + 2a 2 + — + Pok ) k x k y = D yx (k) 



where we have defined a Q = Ai, /3q — A 2 , and 



a, = £ A (A') T (K) , 



(A7) 
(A8) 
(A9) 

(A10) 



« 2 - 2^ t To W Ai 0*0 + Ti 0*0 A ° 

K#0 



(All) 



a 3 = 22 t A ° (*) T2 ( X ) + A2 (*0 T ° K l 



(A12) 



a 4 = t A ° W T 3 W + A i ( K ) T i (*Q + A 3 (K) To (AT)] K 4 X 

K^O 

With the dynamical matrix written as 

D at p (k) = n^ 1 (A + /Lt) k a kjj + nj 1 ^ik 2 8 a ^, 



(A13) 



(A14) 
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the Lame coefficients are thus given by 

A= I — I — +a Q + a 1 + 2a 2 -a 3 + — + (3 k \ , (A15) 

2Trn 2 e 2 \ ( < , > 

-^-](a 3 + T ). (A16) 

The correction f3 k to the Lame coefficient A (and so to the bulk modulus) comes from the K = term in the 
summation over K in Eq. (|A1|I i.e. it comes from the long-range part of the Coulomb interaction. 
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